The NGRDI was used to create the new and reclassified rasterlayers.
Filter: dissimilarity
window size: 3x3, 9x9, 15x15, 31x31, 45x45
library(sp)
library(raster)
library(rgdal)
library(glcm)
rasterOptions(maxmemory=1e+08)
a=raster("C:/master-prog/data/aerial/aerial_croped/NGRDI.tif")
diss_31=glcm(a,window = c(31,31), statistics = "dissimilarity")
writeRaster(diss_31, "C:/master-prog/data/aerial/aerial_croped/NGRDI_31", format="GTiff")
diss_3=glcm(a,window = c(3,3), statistics = "dissimilarity")
writeRaster(diss_3, "C:/master-prog/data/aerial/aerial_croped/NGRDI_3", format="GTiff")
diss_9=glcm(a,window = c(9,9), statistics = "dissimilarity")
writeRaster(diss_9, "C:/master-prog/data/aerial/aerial_croped/NGRDI_9", format="GTiff")
diss_15=glcm(a,window = c(15,15), statistics = "dissimilarity")
writeRaster(diss_15, "C:/master-prog/data/aerial/aerial_croped/NGRDI_15", format="GTiff")
diss_45=glcm(a,window = c(45,45), statistics = "dissimilarity")
writeRaster(diss_45, "C:/master-prog/data/aerial/aerial_croped/NGRDI_45", format="GTiff")
Dissimilarity is a distance-based approach on a “square”-matrix to show distribution between point A in dependence to point B and C. It is commonly used in socio-ecological studies to predict distribution in a space between two species.
library(sp)
library(raster)
library(rgdal)
## rgdal: version: 1.2-15, (SVN revision 691)
## Geospatial Data Abstraction Library extensions to R successfully loaded
## Loaded GDAL runtime: GDAL 2.2.0, released 2017/04/28
## Path to GDAL shared files: C:/Users/Admin/Documents/R/win-library/3.4/rgdal/gdal
## GDAL binary built with GEOS: TRUE
## Loaded PROJ.4 runtime: Rel. 4.9.3, 15 August 2016, [PJ_VERSION: 493]
## Path to PROJ.4 shared files: C:/Users/Admin/Documents/R/win-library/3.4/rgdal/proj
## Linking to sp version: 1.2-5
plot(raster("C:/master-prog/data/aerial/aerial_croped/NGRDI_3.tif"))
library(sp)
library(raster)
library(rgdal)
plot(raster("C:/master-prog/data/aerial/aerial_croped/NGRDI_15.tif"))
library(sp)
library(raster)
library(rgdal)
plot(raster("C:/master-prog/data/aerial/aerial_croped/NGRDI_45.tif"))
Best window size is 15x15 since lower window size (<10) creates to much noise in the raster and bigger than this size (>20) merges to much unknow pixel into one value.